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Abstract 

Context. The QUBIC collaboration is building a bolometric interferometer dedicated to the detection of B-mode polarization fluctu- 
ations in the Cosmic Microwave Background. 

Aims. We introduce a self-calibration procedure related to those used in radio-interferometry to control a large range of instrumental 
systematic errors in polarization-sensitive instruments. 

Methods. This procedure takes advantage of the fact that in the absence of systematic effects, measurements on redundant baselines 
should exactly match each other. For a given systematic error model, measuring each baseline independently therefore allows to write 
a system of nonlinear equations whose unknowns are the systematic error model parameters (gains and couplings of Jones matrices 
for instance). 

Results. We give the mathematical basis of the self-calibration. We implement this method numerically in the context of bolometric 
interferometry. We show that, for large enough arrays of horns, the nonlinear system can be solved numerically using a standard non- 
linear least-squares fitting and that the accuracy achievable on systematic effects is only limited by the time spent on the calibration 
mode for each baseline apart from the validity of the systematic error model. 

Key words. Instrumentation: Polarimeters - Instrumentation: Interferometers - Methods: Data Analysis - Cosmology: Cosmic 
Background Radiation - Cosmology: Inflation 



1. Introduction 

The quest for the B-mode of the polarization of the cosmic back- 
ground is one of the scientific priorities of observational cos- 
mology today. Observing this mode appears as the most power- 
ful way to constrain inflation models. However, detecting such 
a weak signal is a real experimental challenge. In addition to a 
high statistical sensitivity (huge number of horns and bolome- 
ters required), future experiments will need excellent quality of 
foreground removal and unprecedented control of instrumental 
effects. 

Currently, most projects are based on the experimental con- 
cept of an imager. A promising alternative technology is bolo- 
metric interferometry. This is the project of the QUBIC instru- 
ment (the QUBIC collaboration 2010). A first module is planned 
to be installed at the Franco-Italian Concordia Station in Dome 
C, Antarctica in 2014. The aim is to combine the advantages of 
an imager in terms of sensitivity with those of an interferometer 
in terms of control of systematic effects. The statistical sensitiv- 
ity of the QUBIC instrument is comparable to that of an imager 
with the same number of horns covering the same sky fraction. 
The full QUBIC instrument (six modules) will comprise three 
frequencies (97, 150 and 220 GHz) and aims to constrain, at the 
90% confidence level, a tensor-to-scalar ratio of 0.01 with one 
year of data. 

The aim of this article is to introduce a new method, spe- 
cific to bolometric interferometry called self-calibration, and to 
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give an example of application with the QUBIC instrument. This 
method allows to control a wide range of instrumental system- 
atic effects. 

This self-calibration technique is based on the redundancy 
of the receiver array (Wieringa 1991). It uses the fact that in the 
absence of systematic effects, redundant baselines of the inter- 
ferometer should measure exactly the same quantity. For a real 
instrument, these measurements will be different because of sys- 
tematics. These small differences can be used to calibrate param- 
eters that characterize completely the instrument for each chan- 
nel and estimate the instrumental errors. 

In the case of a bolometric interferometer, the square horn 
array will provide a large number of redundant baselines. In this 
way, a bolometric interferometer can be self-calibrated thanks to 
a calibration mode during which it will measure separately the 
rihirih - l)/2 baselines or a fraction of the «/,(«/, - l)/2 base- 
lines with «/, the number of horns observing an external polar- 
ized source. 

This method is inspired by traditional interferometry 
(Pearson & Readhead 1984) where signal phases are often lost 
due to atmospheric turbulence. Standard calibration procedures, 
similar to those used for imaging techniques and based on ob- 
servations of an unresolved source whose flux is assumed to 
be known, exist in radio-interferometry. We emphasize that the 
word self-calibration refers, by opposition, to a procedure in 
which no knowledge of the observed source is required (most of 
the time, the object which is scientifically studied is itself used 
as the calibration source). Most of these self-calibration tech- 
niques are based on the evaluation of so-called closure quantities 
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phases or amplitudes. A set of unknown phases can for instance 
be iteratively reconstructed by forming quantities where they are 
nulled (the product of the three visibilities that can be formed 
with three antennas). The use of redundant baselines for calibra- 
tion is by contrast rather uncommon in radio-interferometry - 
see (Wieringa 1991; Noordam & de Bruyn 1982). This comes 
from the fact that most radio-interferometers have very few re- 
dundant baselines; in order to achieve very high angular resolu- 
tion, it is indeed better to arrange a given number of antennas in 
order to optimize the uv-plane sampling, rather than to maximize 
redundancy. 

A new kind of all-digital radio-interferometer the "omnis- 
cope", dedicated to 21 cm observations, has recently been pro- 
posed by (Tegmark & Zaldarriaga 2009a,b); its concept can be 
summarized by the following five steps: 

1) signals collected by antennas are digitized right after am- 
plification, 

2) a temporal Fast Fourier Transform (FFT) is performed in 
order to split them into frequency subbands, 

3) a spatial FFT is performed on each subband set, 

4) the square modulus of the FFT result is computed, 

5) an inverse spatial FFT is performed to recover the visibil- 
ities. 

There are conceptual similarities with the bolometric inter- 
ferometer concepts, since in this latter instrument, steps 2 and 
3 are performed in an analogical way respectively by the beam 
combiner and the bolometers. Because antennas have to be lo- 
cated on a grid in order for the FFTs to be performed, an om- 
niscope will possess many redundant baselines; (Liu et al. 2010) 
have shown how this allows to self-calibrate the antenna com- 
plex gains. In both cases (standard radio-interferometer or om- 
niscope), the aim is to calibrate the complex antenna gains; there 
is actually a mathematical trick to get a linear system of equa- 
tions from which these gains can be obtained. 

This paper is organized as follows. Section 2 introduces the 
self-calibration method with the Jones matrix formalism and the 
Mueller matrix formalism in the case of radio-interferometry and 
omniscope. We show in section 3 how the procedure can be ap- 
plied to the QUBIC bolometric interferometer and finally, we 
describe the possible self-calibration algorithm and its results. 



2. General principle 

2.1. Instrumental systematics modelisation with Jones 
matrices 



after passing through the receiver E r is 



(2) 



where the Jones matrix J is a 2 x 2 complex matrix. It describes 

how the instrument linearly transforms the two dimensional vec- 

— > 

tor representing the incoming radiation field E into the two di- 
mensional vector of the outgoing field E r . 

For an instrument with several components, the Jones ma- 
trix is the product of the Jones matrices for each component. 
For example, the ideal Jones matrix for an instrument in which 
the incident radiation passing through a rotating half-wave plate 
before propagating through the horns is 



J rot J hwpS rot^ horn 

jcos(cjt) - sin{(jL>t)\l\ 
sin(cjt) cos(ojt) 1 10 -1 



cos((Dt) sin< 
-sin(o)t) cos 



in(cot)\h 0\ 
is(o)t)j\0 lj 



cos{2(jjf) sin(2o>t) 
sin(2toi) -cos(2o>t) 



(3) 



where <x> is the angular velocity of the half-wave plate, J rol 
is the rotation matrix, Jhwp is the ideal Jones matrix of the half- 
wave plate and Shorn is the ideal Jones matrix of one horn. 

After passing through the receiver, ideal orthogonal linear 
detectors measure the power in two components 

51 = + Qcos(Aui) + Usin(4cot)) 

1 

5 2 = -(/ - QcosiAwf) - Usin{Aujt)). 

To model systematic errors within a polarization-sensitive 
interferometer, the Jones matrix can be described by introduc- 
ing diagonal terms: the complex gain parameters g x and g y and 
non-diagonal terms: the complex coupling parameters e x and e y 
associated to the orthogonal polarizations 



1 - gx e x 
e y 1 - gy 



(4) 



In this section, we use the notation proposed by (O'Dea et al. 
2007). With the electric field of an incident radiation at a fre- 
quency <p defined as RlEe^'], and choosing two basis vectors 

e* x and ~e y orthogonal to the direction of propagation k , all the 
statistical information is encoded in the coherence matrix C 



C = 



(\E X \ Z ) (E X E*) 
\{E;E y ) (\E y \ 2 ) 

1 ( I + Q U-iV 



(I + Q U-iV\ 

[u + iv i-qJ' { ' 



Systematic errors arising from the half-wave plate and from 
the square horn array can be modeled by : 

1 . a Jones matrix for the half-wave plate 



2. a Jones matrix for the horn i with 1 < i < n n 



ihornj ~ 



1 



e x ,i 



(6) 



where E x , E y are complex amplitudes of the transverse elec- 
tric field E and I, Q, U and V are the Stokes parameters. 

The propagation of a incident radiation E through a receiver 
can be described by a Jones matrix J such that the electric field 



The electric field E propagated through the half-wave plate 
and the horn ;' becomes 



— J rot J hwp$ rot J horn,i E • 



(7) 
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2.2. Instrumental systematics modelisation with Mueller 
matrix 

The Jones matrix expresses the transformation of the electric 
field in the x and y-directions and the Mueller matrix describes 
how the different polarization states transform. The Jones matrix 
is a 2 x 2 matrix whereas the Mueller matrix is a 4 X 4 matrix. 
So, the 4x4 Mueller matrix can be written as the direct product 
of the 2x2 Jones matrices. 

Let us calculate the tensor product of the outgoing field given 
byEq.Q 

"e;.® ~E r = J£®J*£*. (8) 

In general, the direct product (E ® E*) gives the vector C 
defined as 



These vectors could reduce to a scalar or a vector with 2, 3 or 
4 elements depending on the Stokes parameters the instrument 
is sensitive to. One can define a = 1,2,3,4 the number of Stokes 
parameters the instrument allows to measure. 

The «/,(«;, - l)/2 baselines of the interferometer can be clas- 
sified into n + sets sp of redundant baselines (same length, same 
direction) indexed by (3. In the absence of systematic errors, the 
redundant visibilities should have exactly the same values 



V{i,j)« /h Vy = V^ 



(15) 



For a real instrument however, redundant visibilites Vy will 
not have exactly the same values because of systematic errors 
and statistical (photon) noise, and one can write the following 
system of a x «/,(«/, - l)/2 complex equations 



C 



r 

^ A'V 

c yx 

Cyy 



E x 








® 


e; 



(E X E* X ) 

(E X Ey) 

(E y ED 
{E y El) 



(9) 



Using Eq.([8]l, one can write the transmission of the electric field 

through an instrument described by its Jones matrix using the 
— > 

vector C of this electric field 

C r = JE®T~E* = J®J*C. (10) 

Accordingly Eq.lJTji, the polarization state of this electric 
field can be described by the Stokes vector S defined by 



S = 



C +C 

^xx ~ *~*yy 

C -C 

^xx ^yy 

Cxy "t" Cy X 

i(C x y Cy x ) 



(ID 



where I, Q, U and V are the Stokes parameters. 

One can obtain the expression of the Stokes vector from the 
— > 

vector C defined in Eq.(9) 



S =AC, 



(12) 



where 



10 1 
10 0-1 
110 
i -i 



Substituting Eq.( 12 1 into Eq.( 10 1, it follows that the outgoing 



Stokes vector S r can be written as 



(13) 



where M = A(J ® J*)A 1 is the Mueller matrix which de- 
scribes how the Stokes parameters transform. 

2.3. Polarized measurement equation with Mueller formalism 

A polarization-sensitive interferometer measures the complex 
Stokes visibilities from all baselines defined by the horns i and j 
in an array of receivers 



(V 1 \ 



(14) 



Vh 



Mij.Vp + mj 



(16) 



where n, ; are statistical noise terms and where M/y are a kind 
of complex Mueller matrices which reduce to the identity matrix 
for a perfect instrument. 



My 



(M" M 1Q M ,u M- 



f iv\ 



,Q1 



M Jdu M 'dv 



Ml 1 M? Q MV. u M" v 

.MY 1 mYP MY v MY v 
v IJ IJ I] 'J ■ 



(17) 



The elements of these matrices are not independent and can be 
expressed in terms of the diagonal and the non-diagonal terms 
of the Jones matrix. 

The first order Mueller matrix for a polarization sensitive ex- 
periment is 

1 (A B\ 



Mi. 



2 C D 



(18) 



with 



2 + g xJ + g*yj + g X ,i + gy,i g X J ~ gyj + gx.i ~ , 



yj 
+ . 



2 + g x j + gyj + gx,i + gy 



(19) 



. , + C ; + &xi + e \ i -i(e* r i - e t i - e x i + <? v i)\ 

B= i J I'J ' . l' J • " ' (20) 

» e x.j - e y,j + e x,i ~ <■•:., ~lU\j + <\j ~ (\, ~ 



C = 



lj + e *yJ + e *< + ~ e lj + e yJ ~ ''° + e y* ) 

e*j - e xJ + e yJ ) i(e* x j + e* yJ - e xJ - e yJ )J 



(21) 



D = ( 2+ 8 * xJ + 8 *' J + 8U + 8yjt ' (8 * J ~ 8 * yJ ~ 8x4 + 8y '^ 

V^Sx.j ~ 8yJ + gxj + gyj) 2 + g* xJ + g* y j + g xJ + gy 



(22) 



2.4. Standard radio-interferometry and omniscope 

In the case of instruments that are not designed to measure po- 
larized radiation fields (like most radio-interferometers or om- 
niscopes), the Mueller matrix reduces to its first element (a - 1), 



My-»Jlf£ = (l+ft)(l+*;) 



(23) 



where gi, gj are the complex gains of antennas i and j. A lin- 
ear system can be obtained by taking the complex logarithm of 
Eq.([T6} 

ViJ =v p + d + G* (24) 



3 



M.-A. Bigot-Sazy et al.: Self-calibration: an efficient method to control systematic effects in bolometric interferometry 



and G* 



where v y = log\v uj \, vp = log\v p \, G, = log\gi\ 

log \g* . In order to get an invertible system, one must add an 
overall gain normalization constraint 



n h -\ 



(25) 



i=0 



For a square grid array of «/, > 8 antennas, the problem becomes 
overdetermined since the number of equations is 2 x «/,(«/, - 
l)/2 ~ with n/j(«/i - 1)/2 the number of the measured complex 
visibilities. The number of unknowns is 2n t + 2«/, ~ 6«/, with 
n + the number of complex visibilities of the unique baselines 
and «/, the number of complex gain factors. Things are however 
made subtler by the use of the complex logarithm - see (Liu et 
al. 2010) for details. 



interferometer, the observable is the superposition of the fringes 
formed by the sky electric field passing through a large number 
of back-to-back horns and then focused on the detector plane ar- 
ray. The image on the focal plane of the optical combiner is the 
synthesized image as only specific Fourier modes are selected by 
the receiving horns array. A bolometric interferometer is there- 
fore a synthesized imager whose beam is the synthesized beam 
formed by the array of receiving horns. 

The electric field E'j^hp) collected by the horn i located at 

~~Xi for the polarization r; e {x, y) and measured by a bolometer 
q when all the primary horns are looking at the same radiation 
field E n (h p ) as a function of the observed direction h p is 



(27) 



2.5. Polarization-sensitive interferometer 

In the polarization-sensitive case, on the contrary, Eq.([To*| can- 
not be easily transformed into a linear system. The number of 
unknowns is a x 2n # + 8«/,, while the number of equations is still 
oc n 2 h . In order to get an invertible system, one must add the gain 



normalization constraint given by Eq.(|25[) and another overall 
normalization constraint for the non-diagonal terms of the Jones 
matrix as the leakages are small and randomly distributed 



Tlfi—l 



2>=o. 



(26) 



;=o 



The system (which becomes overdetermined for «/, > 20 if a — 
3) can then be solved using a standard nonlinear least-squares 
method. 



3. Application to the QUBIC bolometric 
interferometer 

3.1. Observables in bolometric interferometry 

In this section, we derive the expression for the power received 
in the focal plane in the case of bolometric interferometry. 

The bolometric interferometer proposed with the QUBIC 
instrument (the QUBIC collaboration 2010) is the millimetric 
equivalent of the first interferometer dedicated to astronomy: the 
Fizeau interferometer. 

The receptors are two arrays of «/, horns: the primary and 
secondary horns back-to-back on a square grid behind the optical 
window of a cryostat. Filters and switches are placed in front and 
between the horn array. The switches will be used only during 
the calibration phase. The polarization of the incoming field is 
modulated using a half-wave plate located before the primary 
horns. This location of the half-wave plate avoids a leakage from 
the Stokes parameter I to the Stokes parameters Q an U if the 
half-wave plate has no inhomogeneities. 

Signals are correlated together using an optical combiner. 
The interference fringe patterns arising from all pairs of horns, 
with a given angle, are focused to a single point on the focal 
plane. Finally, a polarizing grid splits the signal into x and y- 
polarizations, each being focused on a focal plane equipped with 
bolometers. These bolometers measure a linear combination of 
the Stokes parameters modulated by the rotating half- wave plate. 

With an interferometer, the correlation between two re- 
ceivers allows for direct access to the Fourier modes (visibilities) 
of the Stokes parameters I, Q and U. In the case of a bolometric 



where we have introduced the matrices or,„ and Bj(h p ) which 
characterize completely the instrument. 

A matrix a lq is defined for each channel of horns (i) and 
bolometers (q) and includes the geometrical phases induced by 
the beam combiner, the beams of the secondary horns B sec (d q ) 
and the gain g q of the bolometer 



with 



< = 8q j 



. . , la x 

ideal _ f iq 

"'4 " a. 

\ iq 



Xi 



(28) 



B x sec (dq)exp[i2n—d q ]J(v)®(d - d q )dvdd q 
A 



and, 



a>\ q = gq J B^id^expliln^dqjJW&id - d q )dvdd q . 

The detector location is given by the unit vector d„ and A is 
the wavelength of the instrument. The integrations are over the 
surface of each individual bolometer modeled with the top-hat 
like function &(d - d q ) and on bandwidth of the instrument J(v) 
with v the frequency. 

A matrix Bj(h p ) is defined for each channel of pointings (p) 
and horns (i) and includes the primary beam B n prim ffip), the horn 

position Xi and the direction of pointing h p 



Pi (*,) = | - p, h) 



(29) 



with 



and, 



B*(n p ) = B x ■ i (h p )exp(i27Tx i h p ) 



B\(h p ) = B y i (n p )exp(i2nx i h p ). 



The half- wave plate rotates at angular speed d and therefore 
modulates the two orthogonal polarizations such that 



£" A '\ _ lcos(la)t) sin{2u)t) \ IE X 
E' y J ~ I sin(2cof) -cos(2cot) I ' \E y 



(30) 



The electric field collected by one horn i located at x, is 

sKhp) = E^hpmh,,). (3D 
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The electric field collected by one hom i located at x, and reach- 
ing the bolometer q located at d q is 



s'l a (n p ) = E" l {h p )p i {h p )a i q. 



(32) 



The power measured by one polarized bolometer located at 
d q in the focal plane of the beam combiner is then 

S n qp = ^\E'\h p f\B qtS {h p fdn p (33) 

where the synthesized beam B qs (h p ) for the detector q is formed 
by the arrangement of the primary horn array as follows 



Bqjfip) = ^a iq a* jq Bi(h p )B*(n p ). 



(34) 



The synthesized beam depends on the sky direction h p . So, the 
synthesized image is a convolution of the sky and of the electric 
field through the synthesized beam. 



One can rewrite Eq.( 33 i to exhibit the modulation of the po- 
larization induced by the half-wave plate: 

S n qp (t) = S' qp + ^cos(47rcot)Sf D + ^ sin(47Tcot)S^ D (35) 



qp 



'IP 



where e x = 1 for the polarization x, e y — -1 for the polarization 
y and where S qp are the synthesized images on the focal plane 
for each Stokes parameter X = {I, Q, U). The cosine and sine 
coefficients come from the modulation induced by the rotating 
half-wave plate. 

Systematic effects arising at any level of the detection can be 
modeled by associating a Jones matrix to each horn i J/, „ M and 
a Jones matrix of the half-wave plate J/, vvp . 

They can be introduced as defined in Eq.(|7]i. 



3.2. Self-calibration procedure 

During the self-calibration mode, distinct from the ordinary data 
taking mode, the instrument scans a polarized source and mea- 
sures one by one the «/,(«/, - l)/2 synthesized images from all 
baselines or only a fraction of them. In the QUBIC design, this 
can be achieved using switches located between the back-to- 
back horns. The switches are used as shutters that are operated 
independently for all channels and are only required during the 
calibration phase. One can modulate on/off a single pair of horns 
while leaving all the others open in order to access the synthe- 
sized images measured by each pair of horns alone. 

This procedure requires the knowledge of the individual pri- 
mary beams of each horn. The maps of the primary beams can 
be obtained independently through scanning an external unpo- 
larized source. 

By repeating this with all baselines, all bolometers and all 
directions of pointing, one can construct a system of equations 
whose unknowns are 

1. the complex coefficients cP iq , defined for each hom i, each 
bolometer q and for each polarization 77 which correspond to 
4n/,rc 9 parameters, 

2. the horn location x, (2«/, parameters), 

3. the direction of pointing h p (2n p parameters), 

4. the complex horns systematic effects g x j, g y j, e x j, <? V;i - defined 
for each horn i and for each polarization rj (8«/, parameters), 

5. the complex half-wave plate systematic effects h x , h y , % x , g y 
and for each polarization r\ (4«/, parameters). 



For an instrument with rc/, horns, n q bolometers and for a scan 
of n p pointings, the number of unknows is 



n„ = 4 x n h x n q + 2 x n h + 2 x n p + 8 x n h + ! 



(36) 



During the self-calibration procedure, the number of con- 
straints is given by the measurements, i.e. the synthesized im- 
ages. One has the «/,(«/, - l)/2 measured synthesized images 
for each bolometer n q , each pointing n p , each Stokes parame 



ter and each focal plane. With Eq.(25 1 and Eq.po}, one can add 
four complex equations : the sum of the diagonal and the non- 
diagonal terms of the Jones matrices. The number of constraints 
is given by 



6 x n h x (n h - l)/2 xn q xn p + '• 



(37) 



The problem becomes easily overdetermined: for an instru- 
ment with «/, = 9, n q — 4 and n p = 10, the number of constraints 
is 8648 and the number of unknows is 262. It can be solved with 
a least squares algorithm. 

The first module of the QUBIC instrument will consist in 
400 primary horns and two 1024 element bolometer arrays. The 
number of constraints could be reduced if the self-calibration is 
not performed on the «/,(«* - l)/2 baselines but on a fraction of 
baselines. This will be shown in the following. 

However, closing all the switches except two would actu- 
ally change dramatically the thermal load on the cryostat, which 
could affect the bolometric measurements. Fortunately, there is 
a trick explained in Appendix A that allows to indirectly mea- 
sure Si- {d q ) while minimaly changing the thermal load. One 
can show that 



S*. {d q ) = C l (d a ) + C n (L)-2Sl (L) 



(38) 



where S Jj pq {u q ) is the quantity measured by a bolometer q when 
all the switches are open except the i and j, and C 1 ! (a q ), 
CI (a q ) are the powers measured when all the switches are open 
except respectively i or j. Measuring these three terms therefore 
allows to measure (d q ) while keeping the thermal load al- 
most constant. 

However, this also increases the noise. The noise on each 
term is therefore 

ASl pq = ACf pq + AC% q + 4ASl pq 

= 2(n h - \)NET 2 + 4(n h - 2)NET 2 

= (6n h - 10)NET 2 

where NET is the noise equivalent temperature of the bolome- 
ters and T the temperature of the 100% polarized source. 



3.3. Numerical simulation 

We have numerically implemented the method to check if the 
nonlinear system could be solved 

We generate the instrument with a set of ideal parameters 
(horn locations, directions of pointing, primary and secondary 
beams, detector locations...) and a set of parameters randomly 
corrupted by systematic errors (horn location errors, pointing er- 
rors, assymetries of beams, bolometer location errors, diagonal 
and non-diagonal terms of the Jones matrices ...). The widths of 
the random deviation of all corrupted parameters around their 
ideal value are given in Table 1 . 

In order to get a solvable system, one must add some normal- 
ization constraints for the coefficients a n iq and y6j(« p ) which do 
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not change the modelling of systematic errors. They mean that 
the self-calibration only allows for relative calibration of these 
parameters. One can add: 



L. 



an absolute calibration of the global gain of the instrument, 
a" = 1 

a convention on the phase of the a. coefficients, 
Vq, arg(a n 0q ) = 0. 

A rotation of global phase <p q applied to the coefficients a! 
for one bolometer q does not modify Eq.(|27|i and Eq.(|34( 
and therefore the observations, 
a convention on the primary beams, V/,/J?(0)) = 1 
Multiplying the a. coefficients by a term defined as c,e^' and 

dividing the fJ'.(h p ) coefficients by the same term does not 
modify Eq.( 27 1 and Eq.( 34 1 and therefore the observations, 
a convention on the phase of the f3l(h p ) coefficients, 
Vp,arg(J? (h p )) = 0. 

A rotation of global phase (f>(fip) applied to the f3''(h p ) co- 
efficients for one pointing p does not modify Eq.(27i and 
Eq.( 34 1 and the observations. 



We compute the corrupted synthesized images and add gaus- 
sian statistic noise given by 



NETp 
V^V2 



-\/6«/, - 10 



(39) 



where the noise equivalent temperature of the bolometers 
NET is taken to be 3Q0p.K.s*, the temperature of the 100% po- 
larized source is T=100K and the time spent on each baseline on 
the calibration mode is t/, = Is. The usual convention is to give 
the NET for unpolarized detectors but it is convenient in our case 
to use quantities with polarization. In this case, the NET is given 
by NET = M ^ L - 

We solve the non-linear system with a standard non-linear 
least-squares method based on a Levenberg-Marquardt algo- 
rithm. The ideal coefficients (without systematic errors) are used 
as starting guess for the different parameters. 



Error source 


Gaussian level 


Pointing Uncertainty h p 
Horn location error x, 


1 [deg] 


100 [/mi] 


Detector location error d q 


1 [/mi] 


Bolometer gain g q 


0.01 


Primary beam error B' 1 . .(«„) 


0.01 


Secondary beam error If sec (d q ) 


0.01 


Horn systematics g, hi 


0.0001 


Horn systematics e, t , 


0.0001 


Half-wave plate systematics h,, 


0.01 


Half-wave plate systematics 


0.01 



Table 1: Range for systematic errors for each parameter. 



3.4. Results 

In order to be able to have a large number of realizations, we 
run the simulation for an array of 9 primary horns, 9 bolome- 
ters and 10 pointings for 100 realisations. Fig.l shows the result 



of the self-calibration simulation. The six plots are scatter plots 
of ideal vs. real synthesized images in red and of recovered vs. 
real synthesized images in blue for the X and Y focal planes 
and each Stokes parameter I, Q and U. The synthesized images 
are computed with Eq.(33), the ideal synthesized images are the 
synthesized images without systematic effects, the real synthe- 
sized images are the simulated measurements and the recovered 
synthesized images are computed with the output parameters of 
the self-calibration simulation. The red plots show the corrup- 
tion after adding the systematic effects defined in Table 1. The 
blue plots show that the corruption is solved after applying the 
self-calibration method. 

This method allows access to the systematic effects of the 
horns and of the half-wave plate. It also allows to calibrate the 
parameters a! and f3i(h p ) that characterize completely the instru- 
ment for each channel of pointings, horns, bolometers. 

In running the simulation, one can find that the residual error 
on each output parameter will depend on the number of horns, 
bolometers, pointings, baselines per pointings and on the time 
spent measuring each baseline. Adding more horns, pointings, 
bolometers and baselines increases the mathematical constraints 
on a given measurement, it allows to form new baselines and 
adds redundancy on the horn array. Fig. 2 shows that the residual 
diagonal term error of the Jones matrix of the half-wave plate 
improves as the number of baselines per pointing is higher. This 
result was obtained for a simulation with 9 primary horns, 9 
bolometers and 10 pointings for 40 realisations and for a time 
spent on each baseline tb = 1 s. The baselines measured for each 
pointing are chosen randomly. Similar plots are obtained with 
the other parameters defined in Table 2. 

One can put together these variables and define a power law 
which allows to calculate the residual error on each parameter 
defined in Table 1 : 



1 



1 



CT r — CX 



1 



1 



1 



(40) 



with c a constant, k the exponent of the measuring time per 
baseline tb and a, ft, y, O the exponent of the number of horns 
«/,, bolometers n q , pointings n p and the percent of baselines per 
pointing n\, s - 

The values for each index are summarized in Table 2 for two 
different measuring times per baseline f/, = 1 and % = 100s. 

One can observe in Table 2 that the error on the different 
reconstructed parameters is better when the time spent on each 
baseline is longer. Fig. 3 represents the residual half-wave plate 
gain error as a function of the time spent on each baseline during 
the calibration mode. It shows that the limitation of the accuracy 
achieved on the systematic parameters is given by the time spent 
on calibration mode for each baseline tb- 

The first QUBIC module will contain 400 primary horns, or 
79800 baselines; therefore we need to spend about 22h on cali- 
bration in order to measure all the baselines during one second. 
This lapse of time could however be much reduced with a small 
information loss if the self-calibration procedure was not per- 
formed on all baselines. However, the accuracy on the output 
parameters also depends on the number of baselines per mea- 
surement as illustrated Fig. 2. It will be important to determine 
which is the most interesting strategy for the QUBIC instrument. 

Using the law given by Eq.(|40|>, one can extrapolate the re- 
sult given in Table 2 to the residual error for the QUBIC in- 
strument with 400 horns, 2x1024 bolometers and 1000 point- 
ings for two different measuring times per baseline tb — 1 and 
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Figure 1: Results of the self -calibration simulation for the synthesized images S qp Sq P and S qp for the X focal plane at right and Y 
focal plane at left for an instrument with 9 primary horns, 9 bolometers and 10 pointings for a time spent on calibration mode for 
each baseline of t^ — Is and 100 realisations. These plots represent scatter plots of ideal vs. real synthesized image in red and of 
recovered vs. real in blue. The six plots show the advantage of the self- calibration method. 







h = 


Is 






h = 


100 j 




parameters 


a 


P 


7 


o 


a 


P 


r 


<D 


a" 

•i 

h 


0.74 


0.99 


0.75 


0.84 


0.97 


1.55 


1.22 


1.06 


0.46 


0.37 


0.24 


0.73 


0.75 


0.70 


1.02 


1.09 


Xi 


0.58 


0.66 


0.97 


0.82 


0.76 


1.04 


1.28 


1.18 


g,,(x,) 


0.91 


1.06 


0.83 


0.48 


1.04 


1.52 


0.45 


0.75 


e,,(xi) 


0.77 


1.18 


0.36 


0.65 


1.01 


1.16 


0.63 


0.99 


h,, 


0.55 


0.67 


0.12 


0.59 


0.78 


0.85 


0.25 


0.78 


ft 


0.64 


0.58 


0.11 


0.38 


1.11 


0.84 


0.46 


0.80 



Table 2: Results of the self-calibration simulation for each recovered parameter given in the first column for an instrument with 9 
horns, 9 bolometers and 10 pointings. Following the power law given by Eq. \40\ , one can calculate the exponent for each variable: 
the exponent a for the number of horns, the exponent fj for the number of bolometers, the exponent y for the number of pointings 
and O the number of baselines per pointing. This work lias been done for two different measuring times per baseline tj, — Is and 
tj, — 100s and with 40 realisations of the simulation. 



t\, — 100s. The result is given in Table 3 and shows a very signif- 
icant improvement on the level of the residual systematics after 
self-calibration, even for Is. 



3.5. Finding limits 

Accordingly Eq.(35), in an ideal case, in the absence of system- 
atic effects, the powers measured on the X and Y polarized focal 
planes after demodulation of the half-wave plate are 



(T \ 
Tc 
Us J 



-/// 



IB, 



q,s(n p ) 







B, 
















I 


JJ(hp), 



dvdd q dhp 



(41) 



with To, Tc and T$ refer to a constant, cosine and sine terms 
obtained after the demodulation of the half-wave plate and the 
synthesized beam B q s {h p ). The integrations are performed over 
the bandwidth of the instrument, over the surface of the bolome- 
ter and over the sky direction. 

For a real instrument, in the case where the half-wave plate 
is located before the horns, leakages form Q to U and from U to 
Q appear. Using the self-calibration simulation, one can estimate 
the leakage from Q into U and from U into Q by calculating the 
standard deviation of the difference between the ideal and cor- 
rupted parameters cr 1 . , cr Q , , cr u , , err? and cr u ,® 

r r id-corr' id-corr id-corr' id-corr id-corr 

(without self-calibration) and of the difference between the cor- 
rupted and recovered parameters a~ r ' T ' ' 
°~coiT-rec (with self-calibration) 



<-<>rr-/'(;c' ^"corr-rec ® corr-rec ^^1^ 
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parameters 


t b = Qs 


h = U 


t b = 100.5 




O id—corr 


Q~ corr-rec 


ratio 


@~corr— rec 


ratio 


< 
h p 


0.004 


8.48 x 10 " 


47 


1.87 x lO -6 


2140 


0.15 


1.41 x 10" 3 


106 


3.26 x 10~ 5 


4596 




100. x IQ- 6 


5.86 x 10 5 


17 


2.27 x 10~ 8 


4402 




0.0001 


1.36 x 10 6 


73 


1.22 x 10~ 8 


8182 




0.0001 


1.09 x 10- 6 


92 


1.20 x 10~ 8 


8280 


K 


0.01 


1.18 x 10~ 4 


84 


7.27 x 10~ 6 


1375 


ft 


0.01 


1.24 x 10" 4 


80 


5.81 x 10" 6 


1722 



Table 3: Results of the self-calibration simulation for each recovered parameter (given in the first column) for the QUBIC instrument 
with 400 horns, 2 X 1024 bolometers array, 1000 pointings and all baseline measurements. On the second column, one can see the 
value of the standard deviation between the ideal and corrupted parameters (without self-calibration). The third and fifth columns 
give the value of the standard deviation between the corrupted and reconstructed parameters (with self-calibration) for two different 
measuring times per baseline tb — Is and tb — 100s and at total measurement time fixed to be the same for both cases. These values 
are obtained by replacing in Eq.(38) the values of exponent given in Table 2 applied to the design of the QUBIC instrument. The 
benefit obtained after applying the self-calibration method is given in the fourth and sixth columns: this is the ratio between the 
value of the standard deviation of the third column and the second column for tb — Is and of the fifth and the second column for 
t b = 100s. 



diagonal terms - 
power law - 
input error 




* to 4 



10 s : 



10' 




diagonal terms 
power law 
input error 



Percent of baselines per pointing 



10" to 2 10* 10* 

Time [s] spent for each measurement in calibration mode 



Figure 2: Results of the self-calibration simulation for the diag- 
onal terms of the Jones matrix of the half-wave plate and a time 
spent per baseline tb = 1 s. This plot represents the residual error 
on these parameters as a function of the percent of baselines per 
pointing. The red line represents a power law of the shape — 

^bs 

with c a constant, nb s the number of baselines per pointing and 
<t> the index of the power law given in Table 2. This law gives the 
limit of the relative accuracy that can be achieved on the system- 
atic parameters. The green line represents the input error on the 
systematic effects of the half-wave plate in the simulation given 
in Table 1. 



Figure 3: Results of the self-calibration simulation for the diag- 
onal terms of the Jones matrix of the half-wave plate. We show 
the residual error (in blue) on these parameters as a function 
of the measuring time spent on each baseline on the calibration 
mode. The red line represents a power law of the shape ~ 

with NET the noise equivalent temperature of the bolometers, 
T the temperature of the polarized source and tb the measuring 
time spent on each baseline. This law gives the limit of the rela- 
tive accuracy that can be achieved on the systematic parameters 
according to the measuring time spent per baseline. The green 
line represents the input error on the systematic effects of the 
half-wave plate in the simulation given in Table 1. 



In this case without self-calibration, Eq.(39) becomes 



(To) 


• 


T c 









trf. 

id- cor r 



o 





JJ 





tll-Ci 



(W P )\ 

Q(n P ) 

U(hpV 



In this case with self-calibration, Eq.(39) becomes 

^ 



(T \ 
T c 

>Ts , 



cor r— rec 



O UQ 

cr cr 

corr-rec eorr—rec 

rr u rr QU 

u corr-rec corr-rec/ 



Q(n P ) 
U(h p )j 



(43) 



With Eq.(42 1 and Eq.(43 I, one can observe that instrument 
(42) errors and systematic effects induce leakage from Q to U and 
from U to Q but also alter the polarization amplitude. 



8 



M.-A. Bigot-Sazy et al.: Self-calibration: an efficient method to control systematic effects in bolometric interferometry 



It is interesting to focus on the B-mode power spectrum in 
order to estimate the E-B mixing and to constrain the E-mode 
leakage in the B-mode power spectrum. 

— > 

In general, one can define the vector AS which defines the 
errors on Stokes parameters 



10" 



AS 





(A H 


AlQ 


A,u\ 




A QI 


Aqq 


A QU 




[Am 


Auq 


Auuj 



(44) 



where AS = S rec - S corr . 

where S corr is the vector of measured Stokes parameters and 

S rec is the vector of the Stokes parameters obtained with the 
self-calibration method. For the QUBIC instrument, there is no 
leakage of I to Q and U because the half- wave plate is in front of 
the instrument and so, Eq.(44) becomes 



AS 





^QQ ^QU 



(45) 



where Aqq, Aqu, Auq and Auu, the errors between the synthe- 
sized beam without the self-calibration method and the one after 
applying the self-calibration method. 

Following Eq.(45 1, we can define the error matrix on the Q 
and U Stokes parameters 



1 + e p 
-p 1 + e 



(46) 



where the complex term e changes the amplitude of polarization 
and p mixes both the Q and U Stokes parameters. The justifica- 
tion of this matrix is given in Appendix B. 

To go further, one can estimate the leakage from E to B-mode 
and give a constraint on the value of the tensor-to-scalar ratio r. 

An error in diagonal terms Aqq and Auu will affect the am- 
plitude of the E and B-mode power spectrum. The non diagonal 
terms Aqu and Auq results in a leakage from the E to B-mode 
power spectrum (or the B to E-mode power spectrum). 

To have a constraint on the B-mode, as the E-mode is far 
above that of the B-mode in amplitude, one can use the equation 
derived by (C. Rosset et al. 2010) 



ACl BB = Cl BB ' meas - Cl BB 



2eCl BB + p 2 Cl EE 



(47) 



where ACl BB is the global error on the B-mode power spectrum 
Cl BB and Cl EE is the E-mode power spectrum. One can refer to 
Appendix B for explicit details of this equation. It shows that the 
uncertainty on parameter p must be lower than 0.5% to have a 
leakage from E to B-mode lower than 10% for a tensor-to-scalar 
ratio ofr=0.01 for / <100. 



From Eq.( 47 1, one can estimate the leakage from the E-mode 
into the B-mode power spectrum given by the term p 2 Cl EE . Fig. 4 
represents the error on the B-mode power spectrum as a function 
of the multipoles. The colored solid lines represent the leakage 
after applying the self-calibration method for different measur- 
ing times per baseline from tj, = ls to 10000s. The dashed line 
represents the error on the B-mode power spectrum without self- 
calibration method. The black lines are the initial B-mode spec- 
trums for r = 10 _1 , r = 10~ 2 and r = 10~ 3 . 

The leakage from the E-mode into the B-mode is therefore 
significantly reduced by applying the self-calibration procedure, 
even with a modest Is per baseline (corresponding to a full day 
dedicated to self-calibration). The leakage can be further re- 
duced by spending more time on self-calibration. 



10 s 



10" 



r=10' 



r=10' 3 




10s ■ 
100s 



multipoles 

Figure 4: AC; due to leakage from E-mode for different times 
on measurements per baseline % -Is, 10s, 100s obtained with 
the self-calibration method (solid lines) and without the self- 
calibration method (dashed blue line) for the QUBIC instru- 
ment. The black lines are the primordial B-mode spectrums for 
r = 10" 1 , r = IQ- 2 and r = 10~ 3 . 



4. Conclusion 

Bolometric interferometry differs considerably from standard ra- 
dio interferometry in the sense that its primary goal is not to 
reach a good angular resolution but to achieve high statistical 
sensitivity and good control of systematic effects. In this per- 
spective, redundancy turns out to be the crucial property to fulfill 
these two objectives, as shown in (Charlassier et al. 2009) and in 
the present article. 

In this paper, we have shown that with a polarized calibration 
source and the use of the successive observation of this source 
with all pairs of horns of the interferometer (self-calibration), 
one can have low and controllable instrumental systematic ef- 
fects. Redundant baselines should give the same signal if they 
are free from systematics. By modeling the instrument system- 
atics with a set of parameters (Jones matrices, location of the 
horns, beams), one can use the measurements of the different 
baselines to solve a non-linear system that allows to determine 
the systematic effects parameters with an accuracy that, besides 
the correctness of the modelling, is only limited by the photon 
noise, and hence by the time spent on self-calibration. The more 
horns and bolometers in the array, the more efficient the self- 
calibration procedure. 

The resolution of the system is CPU-intensive for large bolo- 
metric interferometers and should be implemented on massively 
parallel computers in the future. Using simulations with various 
horns and bolometer arrays of moderate sizes, we have obtained 
a scaling law that allows us to extrapolate the accuracy of self- 
calibration to the QUBIC instrument with 400 horns, 2x1024 
bolometer arrays and 1000 pointing directions towards the cali- 
bration source. We find that with a few seconds per baseline (cor- 
responding to a few days spent on self-calibration), the knowl- 
edge of the instrumental systematic effects parameters can be 
improved by at least two orders of magnitude, allowing to min- 
imize the leakage from E into B polarization down to a tolera- 
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ble level. This can be improved by spending more time on self- 
calibration. 

The idea of developing bolometric interferometry was mo- 
tivated by bringing together the imager exquisite sensitivity al- 
lowed by bolometer arrays and the ability to handle instrumental 
systematic effects allowed by interferometers. Bolometric inter- 
ferometers have been shown to have a similar sensitivity to im- 
agers (Hamilton 2008 ; The QUBIC collaboration 2010) while 
we show in the present article that the self calibration allows to 
achieve an excellent handling of systematic effects which has no 
equivalent with an imager. 
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Appendix A: Measuring the bolometer power for 
two opened horns is equivalent to measuring 
the bolometer power when all horns are open 
except the horn i and j 

The power collected by the bolometer q for the opened horns i 
and j without polarization is 

S ij = C, + Cj + 2Re(a lq a) q ti iv l3* q ) (A. 1) 

which can be written as: 



S ij = \p i \ 2 + \p j \ 2 + 2Re{p i p)). 



(A.2) 



The total power measured for all the baselines can be ex- 
pressed as 

s t ot = i z p^ 2 = Z lp,n]2 + Z 2Re (p^- (A -v 

m m mH 



The power measured by a bolometer q for all horns opened ex- 
cept the horn i is 

c-i = I Z pj\ 2 = s <°< ~ Ci ~ Z 2Re (PiPD- < A - 4 ) 



The power measured by a bolometer q for all horns opened ex- 
cept the horn j is 

C-i = iZ^'! 2 = Stot ~ C J~Yj 2Re( -PjP*^ (A -^ 

The power measured by a bolometer q for all baselines opened 
except the baseline formed by the horns i and j is 

S _ y = S tot - d -Cj-J] IMPiPl) ~ Z 2Re (PjPd + ZReiPip]) 

k+i k+ j 

(A.6) 

with 

2Re(p iP *) = S tot +S-ij - S tot + C, + 2Re( P iP* k ) - S t0 ,+ 

k+i 

Cj + ^Reipjpl). 

k*i 

Finally, one can find 

2Re(pip*) = S,„, + S ^ -S tot - C-i - C-j. (A.7) 



So : 



S-ij = C-i + C-j + 2Re( P ip)). 



(A.8) 



Appendix B: Error on E and B-modes power 
spectrum 

One can define the Stokes parameters in spin-2 spherical har- 
monics base 

Q(h) + iU(n) = ^ a ±2 i m±2 Yi m (n) (B. 1) 

lm 

where Q and U are defined at each direction h. 

It is convenient to introduce the linear combinations 



E _ a 2Im + a -2lm 
a lm ~ o 



B - a 2lm - a-2lm 

a lm = i j ' 



One can define the two scalar fields 

CO / 



E(n) = Z Z a fm Y ^ 



1=0 m=-i 
I 



B ^ = Z Z aB lm Y lm(n). 



(B.2) 
(B.3) 

(B.4) 
(B.5) 



1=0 m=-l 



Using the coefficients af m and af m , one can construct the an- 



gular power spectrum Cf and Cf as 



Cf E =<l at | 2 > 



cT =<i < i 2 > . 



(B.6) 
(B.7) 
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One can express the Q and U Stokes parameters as a function 
f the coefficients af m and af m 

Qm = \j][(ial-alUY lm (h)+(-ial-al l ) +2 Y lm (h)] (B.8) 



l,m 



U(h) = ^ J][i(ia? m - afj- 2 Yi m (n) + (-ia? m - af m ) +2 Y, m (h)]. 

l,m 

(B.9) 

Using Eq.(B.4), (B.5), (B.8) and (B.9), one can express the 
coefficients af m and af m as a function of the Stokes parameters Q 
andU 

($ = /M-(g$)^ (BIO) 
where the integration is taken over the whole sky and 



M, 



aim 



-{ +2 Y lm {h) +-2 Y, m (h)) -i{ +2 Y lm {h) -_ 2 Y lm (h))\ 
i(+2Yi m (h) -_ 2 Ffo,(n) -(+ 2 Yi m (h) +_ 2 / ' 

(B.ll) 

In the case of a bolometric interferometer, globar errors on 
synthesized beam will affect the amplitude of polarization and 
mix the Q and U Stokes parameters. 

To model systematic errors, one can introduce a Jones matrix 
which describes the propagation of radiation through a receiver 



E r = JE = 



gx 



1 



■ gy 



(B.12) 



where the gain g v and the leakage e v are complex values. 
From this relation, one can construct the Mueller matrix M 
which tells us how the Stokes vector S transforms 



In this case, Eq.(B.lO) becomes 



7m 



= M alm M s 



Q(n) 
U(h) 



(B.16) 



where the coefficients af^" eas and a B ^ eas include systematic ef- 
fects. 

In terms of power spectra, an error in polarization amplitude 
will affect the amplitude of the the E and B -modes power spectra 
and an error which mixes the Q and U Stokes parameters leads 
to a leakage from the E to B-mode (and from the B to E-mode). 

One can easily find that with systematic effects, the coeffi- 
cients af m and af m can be expressed as 



a 



E,meas _ E 



Im 



= a l,n + < Ea lm-P a l,r. 



a 



B,meas _ B 



(B.17) 

*lm = a ?m + ea L + P a r,n (B.18) 

where e is the error in amplitude and p the error of polarization 
leakage. Using Eq.(B.6) and Eq.(B.7), one can obtain 



£jBB,meas 



Cl EE + 2eCl EE +p 2 Cl BB 



Cl BB + 2eCl BB + p 2 Cl EE 



(B.19) 
(B.20) 

where ci EE,meas and ci BB,meas are the E and B-mode power spec- 
tra including systematic effects, Cl EE and Cl BB the input power 
spectra, the complex term e describes the error of amplitude of 
the B-mode power spectrum and the complex term p results in a 
leakage from the E to B-mode power spectrum. 

To focus on the B-mode, one can define the error on Cl BB 
power spectrum 



BB _ (^jBB,meas 



ACF" = CI' 



-Cl BB = 2eCl BB +p z Cl EE . 



BB 



jEE 



(B.21) 



S r = A(J ® J*)A S = MS (B.13) 
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where A= J j j q and S r is the outgoing Stokes vec- 
i -i 

tor and M = A(J <g> J*)A _1 is the resulting error matrix. 



We are only interested on the Q and U Stokes parameters. In 
this case, the error matrix M becomes 



M = 



M 2 \ M 22 



(B.14) 



The first order of this matrix M is 

M u = l+gx + g y +8x + 8y 
M 2l =e y + e* -e x -e* x 
M n = e x + e* - e y - e* 

M 22 = 1 +gx + gy+g* x + gy 

One can generalize this error matrix for Jones matrix of any 
component and rewrite it as 



M s = 



1 +e 



(B.15) 



where the complex term e describes the error of the amplitude of 
polarization and the complex term p mixes the Q and U Stokes 
parameters. 
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